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ABSTRACT 


Climate change is the momentous and persisting change in 
the world’s temperature, precipitation, humidity, and other 
climatic variables. This study, therefore estimated the 
impact of climatic variations on tomato productivity across 
agro ecological zones of Khyber Pakhtunkhwa, Pakistan. 
Panel data for 28 years (1991-2018) across the six districts 
of the agro ecological was used due to availability of data 
on tomato productivity and climatic variables. Yield of 
tomato, area, maximum temperature and rainfall were 
included in the final estimated model. The results indicate 
that the average maximum temperature and average 
maximum temperature square have a significant impact on 
tomato yield. Average maximum temperature has positive 
coefficient while the average maximum temperature square 
has a negative coefficient. This demonstrates that, at first, 
the tomato yield increases as the temperature rises. It 
reaches the maximum at the critical temperature (34.95°C) 
but shows a decline once the temperature rises from the 
critical value. 
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Introduction 


Climate change is a serious global issue faced by the mankind (Auffhammer, 2018). 
It is the momentous and persisting change in the world’s temperature, precipitation, 
humidity, and other climatic variables (Birtahl et al 2021). Studies have identified 
several human activities, industrial waste, vehicles’ emissions, and overall global 
pollution, among others, as driving forces behind global climate change (Lough, 
2008). The world’s economy depends on three major sectors; agricultural, industrial, 
and services. Out of these three sectors, agriculture is the most vulnerable to climate 
change (Parry, 2019) because of its strong dependence on weather patterns. Climatic 
variations in the globe have a significant impact on this sector of the world economy. 
Crop yields, net farm revenue, and agricultural land productivity have all been found to 
be declining in various parts of the world. Global institutions on climate change aim at 
conducting research on climate change and ensuring a sustainable global environment. 
The IPCC being a responsible institution for climate change has forecasted that the 
global temperature will rise to 3°C by the beginning of the next century (IPCC, 2007). 
According to reports from the IPCC climate change has caused damage on a global 
scale. It is responsible for an increase in cyclone activity, an increase in areas facing 
drought, and an increase in heat waves (Bouwer, 2011). From 2007-2016, global 
agricultural activities have polluted the earth by emitting 13% of CO2 and more than 
80% of nitrous oxide. These emissions rise to about 37% if the emissions from pre- and 
post-production activities are added to them (IPCC, 2019). 


Pakistan ranked the 5th most sensitive country to global climatic variation (Abubakar, 
2020). Reasons behind this include lack of policies, less awareness, and increased 
population (Fahad et al., 2020). Climate change is observed impacting agriculture as 
temperature and precipitation are inputs for agricultural production (Feres, Reis, and 
Speranza, 2008) It is more serious issue for farmers of developing countries (Seo et al., 
2009). Being an agricultural country climate change impact on the agriculture of Pakistan 
is a serious challenge for its people and economy. Climate variability has affected the 
production of major crops in Pakistan. During 2018-19, crops in Pakistan showed a decline 
in growth by 4.43%. Production of sugarcane, the major cash crop of Pakistan showed a 
decline of 19.4%. Whereas, production of cotton and rice showed a decline of 17.5% 
and 3.3%, respectively [Government of Pakistan (GOP), 2019]. This shows that climate 
change harms the productivity of major crops in Pakistan. The impact of climate change 
can be reduced by using adaption and mitigation strategies (Chen and Gong, 2021). 


Tomato (Solanum Lycopersicum) is one of the climate-sensitive vegetables consumed 
by people all around the world. The study of the climatic variations’ impact on tomato 
productivity in Pakistan shows that climate change has a significant role in declining 
the yield and production of this vegetable. Due to this, tomato growing farmers in 
Pakistan usually prefer subsistence farming. The use of pesticides for growing this 
crop contributes to a significant share in the emission of gases. These gases boost 
climate change (Ozkan et al. 2011). The yield of tomatoes in the world was recorded 
at 38,272.40 kg/ha in 2018. According to Food and Agriculture Organization USA is 
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the top yielder with estimated yield of 96,807.90 kgha-1, while in Pakistan the yield is 
9,441.2 kgha-1 which shows that Pakistan has still the potential to increase its tomato 
yield per hectare. This higher difference in the yield is because of several factors like 
choice of farmers for crop and acreage allocation, output prices, rain fall, soil type and 
temperature (Birthal et al., 2021) 


Therise in temperature due to climate change decreases the production of vegetables (Hipro 
and Gebeyehu, 2019). Changes in climate scenarios also have economic consequences 
for consumers, such as an increase in commodity prices and a decrease in the utility of 
consuming commodities. The price rise is strongly felt in the case of vegetables and fruits 
due to their strong dependence on climatic conditions. Among many other reason climate 
change is considered as the major one that affect the quality of tomato. The other effects 
of climate change on tomatoes are; tip burn and reduced fruit set. Since 2010, the yield of 
tomatoes in Pakistan has been declining. In 2010, a 9% decline in the yield of tomatoes 
was recorded. In the preceding years, yield kept on showing a decline and failed to recover 
(GOP, 2018). The above-mentioned decline in tomato yield could be due to many possible 
reasons but this research attempts to determine the role of climate change in this regard. 
The decrease in yield is causing an imbalance in the supply and demand of tomatoes in 
Pakistan. For food security purposes, there is a need to implement effective policies to 
maintain the yield of tomatoes in the country (Ahmad and Farooq, 2010). The equilibrium 
in the tomato market would ensure future food security and would bring a fall in the 
import of tomatoes. This study aims to examine the impact of climate variations on tomato 
productivity across agro-ecological zones of Khyber Pakhtunkhwa province of Pakistan. 


Materials and Methods 


Universe of the Study and Data 


The universe of the study of this research was Khyber Pakhtunkhwa (KPK), Pakistan. 
KPK is the third largest province in the country in terms of population. Most of the people 
living in KPK are associated with agriculture for earning their livelihood. The province 
has an influential geographical position that makes it a hub of trade and agriculture. It has 
Gilgit-Baltistan in its north, Afghanistan in the west, Kashmir in the east, and Punjab in 
the south. KP is well-known for the production of many agricultural commodities because 
of different climatic conditions across the province from northern to southern areas. 
The province is divided into 4 agro-ecological zones by the Environmental Protection 
Agency of this province [Environmental Protection Agency (EPA), 2016]. These zones 
are; Northern (A), Eastern (B), Central (C), and Southern (D). Districts from zone A, C, 
and D are selected based on tomato productivity and availability of data. Panel data for 28 
years (1991-2018) was used for the study. Data on climatic variables; average maximum 
temperature, average maximum temperature square, average rainfall and average rainfall 
square was gathered from the Regional Meteorological Department (RMD) Peshawar. 
The data on production, area, and yield of tomatoes in selected districts was taken from 
Crop Reporting Services (CRS) government of Khyber Pakhtunkhwa Pakistan. 
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Conceptual framework 


Several research studies have analyzed the relationship of climate change with crop 
choices, acreage response and its influence on crop comparative advantage (Wang et 
al., 2010; caho and McCarl, 2017; Birthal et al., 2021). To analyze the impact of climate 
change on the productivity of any cereal crop, fruit, or vegetable, a researcher can use 
cross-sectional, time series, or panel data as evident from the literature. The selection 
of the type of data depends on the objectives of the research (Guiteras, 2009). In the 
current research panel data was used to determine the impact of climate variability on 
tomato productivity. The general model for panel data can be presented as; 


Vir = @+ Brie +Eie (1) 


Where in the model Y is the dependent model, X represents the various variables, 1 and 
t represent cross section and time whereas @,8 and B are used for intercept, coefficient 
and error term respectively. Much of the confusion about method of analyzing panel 
data arises due to the fact that different discipline tend to produce solutions according 
to their unique features. This resulted in an astonishing series of notational orthodox, 
terminological variant and various software used. Depending on the researcher’s 
background various models are used for the analyses of Panel data. In literature the 
most widely used model for the panel data are Fixed effect model and Random effect 
model (Niekerk et al., 2022 ). A detailed discussion of the advantages of the both model 
can be found in (Gujarati, D.N., and D.C. Porter. 2009). Fixed-effects model is also 
called an unobserved effect model. It can be presented as; 


yit = (a+ pi) + BXit+Vit (2) 


(a+ pi) in the model shows that the intercept is time-invariant. It means that ,1i will only 

change for its district. The purpose of introducing intercepts in this model is to control 
the time-invariant features (Torres-Reyna, 2007). There are some limitations while 
using the fixed-effect model. These are; the intercepts created for every section that 
would require a degree of freedom and creating dummies increases the possibility of 
strong multicollinearity (Gujarati and Porter, 2009). However, the random effect model 
doesn’t introduce fixed constants for units or sections and considers random parameters 
as intercepts of sections. The Random-effects model forms with the assumption of non- 
correlation of error term with independent variables (Wooldridge, 2013). One simple 
way of understanding the random-effects model is to consider it as a regression model 
with a random intercept or constant (Elhorst, 2014). The generalized form of random 
effects model can be written as; 


yit=a+t+ BXit+(mit+Vit) (3) 
In this general form of the model, the i is time-variant. It means that there will be no 
separate intercepts for districts (Bell, and Jones, 2015). In this paper based on the result 


of the Durbin Wu Hauman test we will decide to use the fixed effect or the random 
effect model (Gujarati and Porter, 2009). 
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Empirical Model 


Several studies have examined the effect of climate change on agriculture. In order to 
find out the impact of climatic variations on tomato productivity across agro ecological 
zones econometric analysis was performed following the existing literature (Bouwer, 
L.M. 2011; Cho and McCarl, 2017; Auffhammer, M., 2018; Hipro and Gebeyehu. 
2019; Chen, and Gong. 2021). Kuamr and singh (2014) estimated that due to rise 
in temperature of about 2.3 OC to 4.5 OC in 2070 to 2099 the food crop grown will 
be declined by 4 percent to 12 in South Asia and Sub-Saharan Africa. The impact 
of climatic and non-climatic on agriculture production have been assessed by many 
empirical studies in the world (Afrin et al., 2017; Chandio et al., 2020; Chao et al., 
2014; Omoregie et al., 2018; Van et al., 2018 and Agbodi et al., 2019). Sarkal et al., 
(2014) conducted a study in Bangladesh and revealed that maximum and minimum 
temperature affects the productivity of agriculture crops. Chandio et al., (2021) also 
pointed out that along with temperature, rainfall, flood, solar radiation and drought 
have an adverse effect on the agricultural productivity. Variation in precipitation and 
temperature adversely affect resources of water and land, which heavily affect the 
agriculture productivity negatively (Mahmood et al., 20122). Ahmad et al., (2020), 
ahsan et al (2020). and Pickson et al., (2020) proposed in their study that rainfall and 
temperature are suitable proxies for climate change. Due to the availability of the data 
on rain fall and temperature the proposed model for estimation is provided as: 


Lny = B0 + B1lnarea + £62 lnmaxtemp + £3 Inmaxtemp2 + 


B4 Inrainfall + B5 Inrainfall2 + Uit 2 


Where Iny is the dependent variable which show tomato yield in kgha-1, Inarea is 
the area under tomato production in ha, B’s are the expected coefficients that need to 
be estimated, Maxtemp and Maxtemp2 represents the average maximum and average 
maximum temperature square, Rainfall and Rainfall2 shows the average rainfall and 
rainfall square, In is the natural log, U is the error term while I and t represents the cross 
section and time period for the study. 


Model Diagnostics 


Several test are used to test which model fits the data well. In the case of the fixed 
effects model, there is a need for several model diagnostic tests as there are chances 
of several issues like; the problem of cross-sectional dependence, heteroscedasticity, 
and time effect autocorrelation (Bakirtas and Akpolat, 2018) in the fixed-effects model. 
To check these issues, different tests were employed. Pesaran test was employed 
to test the first issue. This test was conducted using different statistical software 
(De Hoyos and Sarafidis, 2006). If the p-value of the test is significant there exists 
a problem of cross-sectional dependence. The second problem in the fixed-effects 
model could be heteroskedasticity. This problem was identified by the Wald test of 
group-wise heteroskedasticity. The significant p-value in results shows that group-wise 
heteroskedasticity exists. The third possible issue in the fixed-effects model is of time 
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effect autocorrelation. This was checked by employing a test called Wooldridge test of 
serial autocorrelation. The significant p-value shows the presence of first-order serial 
correlation in panel data (Drukker, 2003). 


With the development in panel data analysis has pointed out the need for cointegration 
and stationarity tests in panel data. The stationarity tests of panel data are more advanced 
than the tests of time series (Bouwer, 2011). The reason is the heterogeneity factor in 
panel data. Some of the simple panel unit-root tests are given. The Im, Pesaran, and 
Shin (IPS) and Breitung tests are two common methods for determining panel data 
stationarity. Due to the problems of cross-sectional dependence, heteroskedasticity, and 
serial correlation in the data set, the panel corrected standard errors (PCSEs) model 
was used to analyze final results. This is an advanced form of fixed effects model and 
estimates results considering the above-mentioned issues in the data. 


Results and Discussion 


In this chapter summary statistics of variables used in the model, results of the analysis 
derived through utilizing the panel corrected standard errors model (PCSE), and critical 
temperature and its impact on tomato yield are presented. 


Impact of Climatic Variations on Tomato Productivity 


Summary Statistics of Variables 


The summary statistics of the variables used in the model is provided in table 1. The 
total observations are 168 i.e., N=6, T=28 and N*T=168. The mean log of yield was 
9.00 kg/ha with a standard deviation of 0.55. The Second variable log of area was in the 
range of 3.40-8.53 ha with a mean of 5.68 and a standard deviation of 1.42. 


Table 1. Summary statistics of variables 


[Variables [Obs [Mean __[Std.Dev. [Min [Max 


9.00 0.55 6.00 9.59 
Iny (Yit) 168 {(9207.70) _|(3645.74)_ | (3272.73) _ | (14620.69) 
5.68 1.42 8.53 
Inarea (Areait) 1681 (363.56) _|(1452.74) [2-408 9-9) | (5051.00) 
3.55 0.07 3.66 
Inmaxtemp (Maxtemp) 168 (34.81) (2.21) 3.35 (28.58) (38.93) 


12.59 0.46 11.24 13.41 
2 
ponaetony: (Maxey?) 168 }(4216.56)  |(150.10) | (816.82) (1515.16) 
3.90 0.63 2.26 5.29 
Inrainfall (Avrainfall) 168 (59.60) (35.32) (9.63) (198.30) 
. 15.64 4.84 27.98 
Inrainfall2 (Avrainfall2) 168 (4792.17) (5827.25) 5.13 (92.74) (3932289) 


Source: Estimated from data, 1991-2018. 


The mean of Inmaxtemp was observed to be 3.55 with a standard deviation of 0.07. 
Log of maxtemp?’ is the fourth variable of the model. Its mean was observed to be 
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12.59 i.e., in the range of 11.24-13.41 and its standard deviation was 0.46. Inrainfall 
is the model’s fifth variable. Its mean is 3.90 with a standard deviation of 0.63. The 
range of mean is 2.26-5.29. The Log of rainfall square is the last variable of the 
model. It has a mean of 15.64 with a standard deviation of 4.84. 


Panel Unit Root Tests 


Panel unit root tests were performed to ensure data stationarity. Stationarity was checked 
for all the six variables used in the model. IPS and Breitung panel unit root tests were 
utilized for this purpose. Stationarity in yield was tested using IPS and Breitung test. 
The results showed that the yield was non-stationary when tested with trend. While it 
became stationary when intercept was added to the trend. It means that the significant 
p-values were obtained when tests included intercept of the yield along with its trend. 


Table 2. Panel Unit Root Tests 


With Trend and Intercept 
PS 


Level 


Temperature : 
° Breitung -5.1154 0.0000*** -6.3217 
CC) 
Average Max. -6.0059 0.0000*** -6.4878 
Temperature ae 
Square (°C) -5.1156 0.0000 -6.3115 
I -6.7616 0.0000*** -6.8424 
-6.2615 0.0000*** -6.3718 


IPS -6.8782 0.0000*** -6.9517 
-6.2545 0.0000*** -6.3597 


Source: Estimated from data, 1991-2018. 
Note: level of significance, ***p<0.01(1%), **p<0.05(5%) 


For checking stationarity in the area i.e., the second variable of the model, both IPS 
and Breitung test were utilized. The result of both tests showed that the area is non- 
stationary when tested with trend only. However, adding intercept with trend made the 
variable stationary. The p-value of IPS for trend and intercept is highly significant and 
shows that the variable is stationary. Average maximum temperature is the third variable 
of the model. Both of the tests of panel unit root show that this variable is stationary 
with trend and with trend and intercept. P-values obtained are highly significant in the 
case of this variable. The same is the case with the fourth variable of the model, i.e., 
maximum temperature square. Both the IPS and Breitung tests show that rainfall is 
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stationary. The p-values obtained for trend and trend and intercept are highly significant 
i.e. 0.0000. The values obtained for rainfall square are also significant and show that 
the variable is stationary. 


Test for Cross-Sectional Dependence and Serial Correlation 


The cross-sectional dependence of data was checked using Pesaran’s test of cross- 
sectional dependence having value of 4.000 with Pr value of 0.0001 indicating the 
existence of cross sectional dependence in the data. To check the heteroscedasticity 
problem Wald test was also employed and the results obtained show that [chi? (6) 
= 3174.45] with P value of 0.000 showing highly significance and the presence of 
heteroscedasticity problem in the data. Panel data autocorrelation was also tested using 
the Wooldridge test for autocorrelation. A highly significant p-value (Prob>F = 0.000) 
was obtained and it was observed that autocorrelation exists in data. 


Estimates of PCSE Model for Panel Data (1991-2018) 


In order to decide which model will be suitable for our data set, we used the Hausman 
test to decide between fixed effects model and random effects model. The result 
obtained from this test shows (Prob>chi2 = 0.0000) highly significant p-value and that 
the best fitted model is fixed effects model. For this study Panel corrected standard 
errors (PCSEs) model was used to analyze final results. This is an advanced form of 
fixed effects model and estimates results considering the issues in the data. Table 3 
shows results for variables used in the model. Area is the first variable used in the study. 
Results reveal that a 1% increase in area will have a positive impact on productivity 
as it will rise by 0.12%. The second variable used in the model was average maximum 
temperature. The significant p value shows that the average maximum temperature 
has a significant impact on the productivity of tomato. The Positive coefficient shows 
that this impact is positive i.e., an increase in temperature will the result in increase in 
tomato productivity. 


By looking at the position of temperature in the table, it can be interpreted that the 
average maximum temperature and average maximum temperature square are 
significantly affecting the tomato yield. The average maximum temperature has positive 
coefficient while the average maximum temperature square has a negative coefficient. 
This implies that the tomato yield initially increases as the temperature rises. It reaches 
the maximum at the critical temperature but shows a decline once the temperature 
rises from the critical level. Pefia and Hughes (2007) Lipper et al. (2009), Shakoor et 
al. (2011), Loum and Fogarassy (2015), and Ghalib et a/. (2017) and Hamdullah et al. 
(2020) find out similar finding in their studies. The average rainfall and average rainfall 
square have an insignificant impact on the productivity of tomatoes in districts studied 
for this research. Our results are in line with results estimated by Islam et al. (2009), 
GCISC (2009), and Khan et a/. (2018). The reason behind the insignificant impact of 
rainfall on tomato productivity is the fluctuation observed in the rainfall pattern due to 
climate change. 
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Table 3. Estimates of PCSE model for panel data (1991-2018) 


fas __—‘ianber ofa fe 

(a 
Estimated coefficients Prob>chi2, fos 
Panel-corrected 


Source: Estimated from panel data, 1991-2018. 
Note: level of significance, ***p<0.01(1%), **p<0.05(5%), ns shows non-significant 


Variation in Yield in Response to Change in Temperature 


The critical temperature for tomato productivity was also calculated from the values 
given in table 3 by using the following formula. 


Critical temperature = exp (B,/2* B,) (5) 
= exp((-129.4095)/(2*(-18.20628))) = 34.95 °C 


The critical temperature for the province indicates that tomato yield will be highest in 
districts where the temperature reaches 34.95°C during the vegetable’s kharif cropping 
season. The maximum yield at critical temperature for the province is estimated to 
be 9763.050kh/ha. However, the yield will start declining when the temperature 
increases this critical value. A graph provided below was constructed to illustrate this 
relationship. Based on this graph and estimated critical temperature for the province, 
the critical temperature for every district used in the study was also estimated. 
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Yield in Kgs 


Figure 1. Variation in yield in response to change in Temperature. 
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Figure 2. Variation in yield in response to change in temperature for districts 
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Source: Authors’ estimated from PCSE model for panel data, 1991-2018. 


Variation in Yield in Response to Change in Temperature for Districts 


Figure 2 represents maximum yield points for their respective districts. The maximum 
yield at critical temperature (34.95 ° C) in Dir Lower will be 9245.93 kg/ha. The yield 
is estimated to decline to 9100.52 kg/ha with a rise in temperature to 36 °C. For Swat, 
the maximum yield at critical temperature is estimated to be 12143.47 kg/ha. The graph 
shows that this yield will decline to 11952.49 kg/ha with rise in temperature to 36 °C. The 
maximum yield for Nowshera at 34.95 °C is estimated to be 7416.12 kg/ha. This yield will 
decline to 7299.49 kg/ha with a 1°C rise in temperature. In Charsadda, the maximum yield 
is estimated to be 8871.42 kg/ha. A decline to 8871.42 kg/ha is estimated with 1°C rise in 
temperature. In D.I. Khan this yield is estimated to be 7530.69 kg/ha at critical temperature. 
It is expected to decline by 118.43 kg/ha with a 1°C rise in temperature. While in the Tank 
the maximum yield is expected to be 7590.62 kg/ha. With a 1°C rise in temperature this 
yield will decline by 119.37 kg/ha. All the graphs show that the yield for respective districts 
is highest at critical temperature and starts to decline with rise in the temperature. 


Figure 3. Zone-wise variation in yield in response to Temperature 
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Source: Authors’ estimate from PCSE model for panel data, 1991-2018. 
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Forecasting the Impact of Rising Temperature on Yield of Tomato According to 
Different Climate Change Scenarios 


To study the climate trend in the country, the meteorological department of Pakistan 
and the Global Change Impact Studies Centre (GCISC) have carried out several studies 
(slam et al. 2009). Tables 4 and 5 present a future scenario for temperature rise in 
relation to tomato impact based on their research. Table 5 shows the response of yield 
to rise in temperature more than 34.95 °C 1.e., critical temperature. The table shows that 
one degree rise in temperature will lower the yield of selected districts by 1.57%. This 
means that tomato growers across zones will suffer losses in terms of yield. The table 
also explains the impact of a 2 °C rise in temperature on yield. It is estimated that this 
rise in temperature will lower tomato yield in selected districts by 5.73%. This reveals 
that climate change in the long-run is significantly harmful to tomato yield across agro- 
ecological zones of Khyber Pakhtunkhwa. 


Table 4. Forecasting the impacts of temperature rise on tomato yield in Khyber Pakhtunkhwa 
according to different climate change scenarios 


Climate change scenarios Yield (Kg/ha) Yield (Yoage) 


Temperature increase by 1 °C -153.54 
Temperature increase by 2 °C -559.26 


Source: Authors’ estimate from PCSE model for panel data, 1991-2018. 


Table 5. Zone wise forecasting the impact of rising temperature on yield of tomato 
according to different climate change scenarios 


£2 


Diktan 
Tink 


Source: Authors’ estimate from PCSE model for panel data, 1991-2018. 


Charsadda -141.75 -1.57% 


Lowest and highest maximum temperature for districts 


Table 6 shows the lowest maximum and highest maximum temperature for districts 
used in the study. The highest maximum temperature in Dir Lower shows that tomato 
yield in the district will increase with the rise in temperature as the district’s temperature 
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hasn’t crossed the critical temperature value i.e. 34.95°C. However, temperature rise in 
the second district of this zone i.e. Swat will cause a decline in tomato yield as the 
district’s highest maximum temperature has already crossed the critical value of 34.95 
°C. Nowshera and Charsadda, are already experiencing losses in the yield due to the 
highest maximum temperature of 38 .93 °C. Also, districts taken from zone D have the 
highest maximum temperature more than the critical value calculated in the study. This 
means that in these four districts the tomato yield is declining. 


Table 6. Lowest and highest maximum temperature for districts 


Lowest maximum temperature Highest maximum temperature 
Dir Lower 28.58 °C 32.08 °C 


Swat 31.86 °C 35.42 °C 
Nowshera 34.6 °C 38.93 °C 


Charsadda 34.6 °C 38.93 °C 


Source: Authors’ estimate for panel data, 1991-2018. 


Conclusions and Recommendations 


This study is aimed at estimating the impact of climate change on tomato productivity 
across agro ecological zones of Khyber Pakhtunkhwa (KP), Pakistan. Three agro 
ecological zones of KP i.e. A, C, and D were selected for this study. Districts from each 
zone are taken based on tomato productivity and data availability. The total districts 
are six i.e. Dir Lower and Swat from zone A, Nowshera and Charsadda from zone C, 
and D.I. Khan and Tank from zone D. Panel data was used for studying the impact 
and following variables. Secondary data for both climatic and non-climatic variables 
were used. Data on climatic variables was collected from the Regional Meteorological 
Department (RMD) Peshawar. While data on non-climatic variables was gathered from 
Crop Reporting Services (CRS). Fixed effects model was selected based on the results 
of the Hausman test. The Data set was also tested for contemporaneous correlation, 
heteroskedasticity, serial correlation, and stationarity. Results revealed that cross- 
sectional dependence, Heteroskedasticity, and serial correlation exist in data. Final 
results were estimated using panel corrected standard errors (PCSEs). The average 
maximum temperature and average maximum temperature square have a significant 
impact on tomato yield. The average maximum temperature has a positive coefficient 
while the average maximum temperature square has a negative coefficient. This implies 
that the tomato yield initially increases as the temperature rises. It reaches the maximum 
at the critical temperature but shows a decline once the temperature rises from the 
critical level. The average rainfall and average rainfall square have insignificant impact 
on the productivity of tomato in districts studied for this research. These results are 
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in line with results estimated by Islam et a/. (2009), GCISC (2009), and Khan et al. 
(2018). The reason behind the insignificant impact of rainfall on tomatoes productivity 
is the fluctuation observed in the rainfall pattern due to climate change. 


The Critical temperature for the maximum yield of tomato was calculated to be 
34.95°C. It was estimated that the yield in all districts of the study showed a decline 
with a rise in temperature above critical temperature. Based on the results of this study, 
it is recommended that policymakers should encourage tomato growers in Dir Lower 
and Swat so that yields in these districts can be increased. To reduce temperature 
rises, the government must concentrate and accelerate tree planting in districts such 
as Charsadda, Nowshera, D.I. Khan, and Tank. Also, heat resistant varieties of tomato 
should be developed for farmers of Nowshera, Charsadda, D.I. Khan and Tank to cope 
with the increase in temperature. Extension officers have to provide required guidance 
to farmers producing tomatoes in selected districts regarding climate change and its 
impacts. 


Limitations of the Study 


The study has several limitations. District Mansehra is top produced tomatoes from 
zone B, but this district was omitted due to statistical discrepancies in data on area 
under tomato and tomato production. Therefore, it can be said that the study doesn’t 
cover all four agro-ecological zones of the province. Individual dummies for districts 
were not created because of the high variation in the data. The model used for the study 
used four basic climatic variables i.e., maximum temperature, maximum temperature 
square, rainfall, and rainfall square. Other important climatic variables could also be 
used for study. These include; minimum temperature, humidity, and sunshine. 


Conflict of interests 


The authors declare no conflict of interest. 


References 
1. Abubakar, S. M., (2019, December 4). Pakistan Sth most vulnerable country to 
climate change, reveals Germanwatch report. Dawn. 


2. Agboklou, K. E., & Ozkan, B. (2023). Analysis of Rice Farmers’ Access To 
Agricultural Credit And Profitability In South Togo . Journal of Business 


3. Economics and Finance.12(1), 1-14. DOI: 10.17261/Pressacademia. 2023. 1722. 


4. Auffhammer, M. (2018). Quantifying economic damages from climate change. 
Journal of Economic Perspectives, 32(4), 33-52. 


5. Ahmad, M., & Farooq. U. (2010). The state of food security in Pakistan: Future 
challenges and coping strategies. The Pakistan Development Review, 49(4), 903-923. 


390 http://ea.bg.ac.rs 


10. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


Economics of Agriculture, Year 70, No. 2, 2023, (pp. 377-394), Belgrade 


Afrin, S., Haider, M.Z. and Islam, M.S. (2017). Impact of financial inclusion on 
technical efficiency of paddy farmers in Bangladesh. Agricultural Finance Review, 
77(4), 484-505. https://doi.org/10.1108/AFR-06-2016-0058. 

Bakirtas, T., & Akpolat, A.G. (2018). The relationship between energy consumption, 
urbanization, and economic growth in new emerging-market countries. Energy, 
147.110-121. https://doi.org/10.1016/j.energy.2018.01.011 

Bell, A., & Jones, K. (2015). Explaining fixed effects: Random effects modeling 
of time-series cross-sectional and panel data. Political Science Research and 
Methods. 3(1), 133-153. 

Birthal, P., et al. (2021). Impacts of climatic hazards on agricultural growth in 
India. Climate and Development. 13: 1-14. https://doi.org/10.1080/ 

Bouwer, L.M. (2011). Have disaster losses increased due to anthropogenic climate 
change? Bulletin of the American Meteorological Society, 92(1), 39-46. 


. Chandio, A., Magsi, H., & Ozturk, I. (2020). Examining the effects of climate 


change on rice production: case study of Pakistan. Environmental Science and 
Pollution Research, 27, 7812-7822. doi:10.1007/s11356-019-07486-9 

Kumar, A. and A. Singh (2014). Climate Change and its Impact on Wheat 
Production and Mitigation Through Agroforestry Technologies. /nternational 
Journal on Environmental Sciences. 5, 73-90. 

Chen, C., Zhou, G.-S., & Zhou, L. (2014). Impacts of Climate Change on Rice 
Yield in China From 1961 to 2010 Based on Provincial Data. Journal of Integrative 
Agriculture, 13, 1555-1564. doi:10.1016/S2095-3119(14)608 16-9 

Chen, S., & Gong. B. (2021). Response and adaptation of agriculture to climate 
change: Evidence from China. Journal of Development Economics. 148, 1-17. 
Cho, S. J., & McCarl, B. (2017). Climate change influences on crop mix shifts in 
the United States. Scientific Reports, 7, 40845. doi:10.1038/srep40845. 

De Hoyos., R.E. & Sarafidis, V. (2006). Testing for cross-sectional dependence in 
panel-data models. The stata journal, 6(4), 482-496. 

Drukker., D.M. (2003). Testing for serial correlation in linear panel-data models. 
The stata journal, 3(2), 168-177. 

Elhorst, J.P. (2014). Spatial panel data models. In Spatial econometrics. Springer. 
Berlin, Heidelberg. 37-93. 

Environmental Protection Agency. (2016). Government of Khyber Pakhtunkhwa 
forestry, environment and wildlife department, Khyber Pakhtunkhwa, Pakistan. 
Ahmad, S., Tariq, M., Hussain, T., Abbas, Q., Elham, H., Haider, I, & Li, X. 
(2020). Does Chinese FDI, Climate Change, and CO2 Emissions Stimulate 
Agricultural Productivity? An Empirical Evidence from Pakistan. Sustainability, 
12(18). doi:10.3390/sul12187485 


http://ea.bg.ac.rs 391 


Economics of Agriculture, Year 70, No. 2, 2023, (pp. 377-394), Belgrade 


21. 


22. 


23. 
24. 


25. 


26. 


27. 


28. 


29. 


30. 


31. 


32. 


33. 


34. 


392 


Ahsan, F., Chandio, A., & Fang, W. (2020). Climate change impacts on cereal 
crops production in Pakistan: Evidence from cointegration analysis. /nternational 
Journal of Climate Change Strategies and Management, ahead-of-print. 
doi: 10.1108/IJCCSM-04-2019-0020. 


Fahad, S., T. Inayat., J. Wang., L. Dong., G. Hu., S. Khan., & Khan., A. (2020). 
Farmers’ awareness level and their perceptions of climate change: A case of Khyber 
Pakhtunkhwa province, Pakistan. Land Use Policy, 96,1-8. 


FAO. 2018. FAOSTATS, (2018) and previous years. www.fao.org 


Feres, J., Reis., E & Speranza., J. (2008). Assessing the impact of climate change 
on the Brazilian agricultural sector. /nstituto de Pesquisa Economica Aplicada 
(IPEA), 1-15. 

Ghalib, H., Shah, S., Jan, A., & Ali, G. (2017). Impact of Climate Change on 
Wheat Growers Net Return in Khyber Pakhtunkhwa: A Cross-Sectional 
Ricardian Approach. Sarhad Journal of Agriculture, 33. doi:10.17582/journal. 
sja/2017/33.4.591.597 


GCISC. (2009). Global Change Impact Studies Centre Islamabad, Pakistan. www. 
gcisc.org.pk 

GOP. Pakistan Bureau of Statistics. (2018). Ministry of Finance, Statistics Division, 
Islamabad, Pakistan. 


GOP. Pakistan Economic Survey. (2019). Economic Adviser’s Wing, Finance 
Division, Islamabad, Pakistan. 


Gujarati, D.N., & Porter., D.C. (2009). Basic Econometrics. Fifth Editon. McGraw 
Hill Irwin USA. 


Guiteras, R. (2009). The impact of climate change on Indian agriculture. Manuscript, 
Department of Economics, University of Maryland, College Park, Maryland. 1-54. 


Tareen, H., Ali, S., Khan, A., Shah, S., & Khan, S. (2021). Economic appraisal 
of transformative climate change on potential variations in wellbeing of wheat 
growers across various ecological zones. Environmental Science and Pollution 
Research, 28, 1-15. doi:10.1007/s11356-020-11409-4 


Hirpo, F. (2019). Review on the Effects of Climate Change Variability on 
Horticultural Productivity. International Journal of Environmental Sciences & 
Natural Resources, 17. doi:10.19080/IJESNR.2019.17.555969. 


Islam, S., Rehman, N., Sheikh, M., & Khan, A. (2009). Climate Change Projections 
for Pakistan, Nepal and Bangladesh for SRES A2 and A1B Scenarios using outputs 
of 17 GCMs used in IPCC - AR4. 

IPCC. (2007). Summary for policymakers. Climate change 2007: the physical 
science basis, Working Group I contribution to IPCC fourth assessment report: 
climate change 2007. Geneva International Development Research — Center, 
Working, Paper No. 3. 


http://ea.bg.ac.rs 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


45. 
46. 


47. 


Economics of Agriculture, Year 70, No. 2, 2023, (pp. 377-394), Belgrade 


IPCC. (2019). Climate Change and Land: An IPCC Special Report on climate 
change, desertification, land degradation, sustainable land management, food 
security, and greenhouse gas fluxes in terrestrial ecosystems. Summary for 
Policymakers. 1-43. 


Khan, A., Ali, S., Shah, S.A., & Fayaz, M. (2018). Impact of temperature and 
precipitation on net revenue of maize growers in Khyber Pakhtunkhwa, Pakistan. 
Sarhad Journal of Agriculture. 34(4), 729-739. 


Lippert, C., Krimly, T., & Aurbacher, J. (2009). A Ricardian analysis of the impact 
of climate change on agriculture in Germany. Climatic Change, 97, 593-610. 
doi:10.1007/s10584-009-9652-9. 


Lough., J.M. (2008). 10th anniversary review: a changing climate for coral reefs. 
Journal of Environmental Monitoring, 10(1), 21-29. 


Loum, A. & Fogarassy, C. (2015). The effects of climate change on cereals yield 
of production and food security in Gambia. Applied Studies in Agribusiness and 
Commerce — APSTRACT. 9(4), 83-92. 


Sarker, M., Alam, K., & Gow, J. (2014). Assessing the effects of climate change on 
rice yields: An econometric investigation using Bangladeshi panel data. Economic 
Analysis and Policy, 44. doi:10.1016/j.eap.2014.11.004. 


Mahmood, N., Ahmad, B., Hassan, S., & Khuda, B. (2012). Impact of temperature 
ADN precipitation on rice productivity in rice-wheat cropping system of Punjab 
province. Journal of Animal and Plant Sciences, 22(4), 993-997. 


Van Niekerk, M., Neser, F., van Wyk, J., & Ducrocq, V. (2022). Comparison of 
fixed and random regression models for the analysis of milk production traits 
in South African Holstein dairy cattle under two production systems. Livestock 
Science, 267, 105125. doi:10.1016/j.livsci.2022.105125. 

Omoregie, O. K., Ikpesu, F., & Okpe, A. E. (2018). Credit Supply and Rice 
Output in Nigeria: Empirical Insight from Vector Error Correction Model 
Approach. /nternational Journal of Economics and Financial Issues, 8(5), 68-74. 
Ozkan, B., Ceylan, R., & Kizilay, H. (2011). Energy inputs and crop yield 
relationships in greenhouse winter crop tomato production. Renewable Energy, 
36, 3217-3221. doi:10.1016/j.renene.2011.03.042. 

Parry, M.L. (2019). Climate change and world agriculture. Routledge. 


R de la, P., & Hughes, J. (2007). Improving Vegetable Productivity in a Variable 
and Changing Climate. An Open Access Journal published by ICRISAT. 4(1), 1-22. 
Seo, S., Mendelsohn, R., Dinar, A., & Kurukulasuriya, P. (2009). Adapting to 
Climate Change Mosaically: An Analysis of African Livestock Management by 
Agro-Ecological Zones. The B.E. Journal of Economic Analysis & Policy, 9, 4-4. 
doi: 10.2202/1935-1682.1955. 


http://ea.bg.ac.rs 393 


Economics of Agriculture, Year 70, No. 2, 2023, (pp. 377-394), Belgrade 


48. 


49. 


50. 


ot. 


52. 


53. 


394 


Shakoor, U., Saboor, A., Ali, I., & Mohsin, A. (2011). Impact of climate change on 
agriculture: Empirical evidence from arid region. Pakistan Journal of Agricultural 
Sciences, 48, 327-333. 

Kumar, A., & Singh, A. (2014). Climate change and its impact on wheat production 
and mitigation through agroforestry technologies. International Journal on 
Environmental Sciences, 5, 73-90. 

Torres-Reyna, O. (2007). Panel data analysis fixed and random effects using Stata 
(v.4.2). Data & Statistical Services, Princeton University, 112. 

Oort, P., & Zwart, S. (2017). Impacts of climate change on rice production in Africa 
and causes of simulated yield changes. Global Change Biology, 24. doi:10.1111/ 
gcb.13967. 

Wang, J., Mendelsohn, R., Dinar, A., & Huang, J. (2010). How Chinese farmers 
change crop choice to adapt to climate change. Climate Change Economics (CCE), 
01, 167-185. doi:10.1142/S2010007810000145. 

Wooldridge, J.M. (2013). Introductory Econometrics A Modern Approach. Fifth 
Edition. South-Western Cengage Learning. 


http://ea.bg.ac.rs 


